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We investigate dynamic scaling properties of the two-dimensional gauge glass model for the vortex 
glass phase in superconductors with quenched disorder. From extensive Monte Carlo simulations we 
obtain static and dynamic finite size scaling behavior, where the static simulations use a temperature 
exchange method to ensure convergence at low temperatures. Both static and dynamic scaling 
of Monte Carlo data is consistent with a glass transition at zero temperature, with correlation 
length exponent given by 1/u = 0.36 ± 0.03. We study a dynamic correlation function for the 
superconducting order parameter, as well as the phase slip resistance. From the scaling of these 
two functions, we find evidence for two distinct diverging correlation times at the zero temperature 
glass transition. The longer of these time scales is associated with phase slip fluctuations across 
the system that lead to finite resistance at any finite temperature. The shorter time scale can be 
described by the form t ~ £*, with a dynamic exponent z = 2.7 ± 0.2, and corresponds to local 
phase fluctuations. 

PACS numbers: 74.25.Qt, 64.70.Pf, 75.40.Mg 



I. INTRODUCTION 

Despite intense studies of the superconducting glass 
transition in vortex systems, the precise nature of 
the transition remains controversial. The presence of 
quenched disorder in a strongly type II superconductor 
in an applied magnetic field may create a superconduct- 
ing vortex glass phase with zero linear resistance, where 
the vortices are localized in random positions given by 
the disorder and interactions 1 .. Considerable theoretical 
and simulation efforts have been focused on whether a 
superconducting glass state exists at finite temperature, 
or alternatively, only at zero temperature. In particular 
the situation in two dimensions still remains debated. In 
this paper we present results from Monte Carlo simula- 
tion and finite size scaling analysis for the glass transi- 
tion in two dimensions. We focus mainly on the dynamic 
scaling properties, which have not been studied in detail 
before, using a dynamic scaling ansatz for a zero temper- 
ature transition. From our static and dynamic scaling 
analysis, we find evidence for a zero temperature glass 
transition, and present various predictions for the physi- 
cal properties at the transitions. 

The gauge glass (GG) model is an XY model with a 
random quenched vector potential on the links of the lat- 
tice; see Eq. Q below. The GG model contains impor- 
tant aspects of order parameter symmetry, interactions, 
and quenched disorder, and has been frequently studied 
as a model for a vortex glass transition. In a real glass, 
the system may fall out of equilibrium at the transition, 
but here we focus on the equilibrium critical behavior of 
the model. This enables the study of universal scaling 
laws at the transition, that depend only on the univer- 
sality class. It has not been settled if the GG model rep- 
resents the correct universality class of the vortex glass 
transition, and it is thus motivated to obtain the precise 
critical behavior of the GG model. The GG model is 
not realistic in a few important aspects: disorder is in- 



cluded only in terms of a random vector potential, instead 
of a random pinning energy for the vortex cores, and it 
has zero applied magnetic field. Alternatively, a suitable 
modification of the GG model as a model for the elastic 
energy of a dislocated vortex lattice is discussed in Rcf. 2. 
Another worry is the use of a classical model for a zero 
temperature transition, where quantum effects may be 
important. However, since vortices localize in the glass 
phase, their quantum fluctuations are small and can be 
neglected, and a classical treatment is a reasonable start- 
ing point, also at T = 0. This neglects the possibility of 
quantum tunneling of vortices, which may actually be 
important. 

Considerable simulation efforts have been invested in 
trying to locate the lower critical dimension d c of the 
vortex glass transition, i.e., the limiting dimensionality 
such that a glass phase exists at finite temperature only 
for d > d c . In the 3D case, several simulations^^ show 
that a vortex glass phase does exist at finite temperatures 
in the absence of magnetic screening. In 2D the situa- 
tion is less clear. It has been shown rigorously that there 
can be no vortex-glass order in 2D for T > For a 
finite screening length it is expected that no glass phase 
can exist at T > 7 . Akino and Kosterlitz computed 
the domain wall stiffness exponent 9 at zero tempera- 
ture, and obtained a negative value, 8 — —0.36, indi- 
cating that T c — 8 . Holme and Olsson computed the 
energy barrier for vortex motion across the system and 
argue that T c > OS. Katzgraber and Young recently ob- 
tained T c = and 9 = -l/v = -0.39 ± 0.03 from finite 
temperature simulations 3 '--. In contrast, several other 
studies of static quantities obtained T c « 0.2 11 ' 12 . Dy- 
namic scaling in the 2D GG model has been studied in 
several papers 12 ' 13 ' 14 ' 1 ^^^, with conflicting results. It 
is important to consider dynamic scaling in order to ob- 
tain information about physically measurable quantities 
like the resistance. Previous studies of zero tempera- 
ture scaling of dynamic properties considered nonlinear 
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I — V characteristics^, and in addition the linear resis- 
tance was found to be roughly consistent with a simple 
Arrhenius form 13,14 , R ~ e~ K / T , but no further support 
for zero temperature dynamic scaling has been obtained. 
In contrast, Refs. Il2lll5lllq obtained T c > from vor- 
tex dynamics simulations. Hence the behavior of the GG 
model remain controversial, in particular for the dynamic 
properties. 

In this paper we study critical scaling properties at 
the glass transition for both dynamic and static quanti- 
ties, assuming a zero temperature transition. Our simu- 
lations are motivated in part by new methods that have 
recently been developed. At low temperature the cor- 
relation times in the simulation grow very long and ef- 
ficient convergence acceleration methods become highly 
desirable. The temperature exchange methodiS* 1 ^ has 
proven very useful for such problems, where configura- 
tions are exchanged between different temperatures, in 
order to avoid getting stuck in metastable states. This 
considerably speeds up the convergence of the simulation, 
and thereby enables simulation at larger system sizes and 
lower temperatures. Our scaling results for static quan- 
tities are in good agreement with some of the previous 
results in the literature for the correlation length expo- 
nent and the domain wall energy exponent #&2ii2ii&ii. 
For our dynamics simulations we also get valuable input 
from the exchange method, even though exchange is not 
directly applicable for dynamics, by using typical equilib- 
rium configurations from the exchange runs as input to 
the dynamics runs. This considerably extends the regime 
of temperatures and system sizes accessible for the dy- 
namics simulation. The main new results of our work 
are for the dynamic finite size scaling analysis of the lin- 
ear resistance and of the dynamic correlation function for 
the order parameter at the zero temperature glass transi- 
tion. We obtain correlation times from finite size scaling 
of our dynamic MC data that show novel behavior. We 
find evidence for two distinct diverging correlation times, 
one associated with phase slip fluctuations, and the other 
associated with local phase fluctuations that are indepen- 
dent of the phase slip processes. 

The organization of the paper is as follows: In Sec. ITT1 
we introduce the GG model and describe our simulation 
methods. In Sec. IHII we discuss the quantities calculated 
in the simulation and the finite size scaling methods. In 
Sec. El the static results are presented, and Sec. con- 
tains the results for the dynamic quantities. Section I VII 
contains discussion and conclusions. 



II. MODEL AND MONTE CARLO 
SIMULATION DETAILS 

The gauge glass (GG) mode l 21 ! 1 ' 21 ! 22 is defined by the 
Hamiltonian 

H = -J2^s(9 i -9 j -A ij ) (1) 

<i,J> 



where Oi is the phase of the superconducting order pa- 
rameter on site i of a square lattice, and Y^Uj) denotes 
summation over all nearest-neighbor pairs in the lattice. 
The vector potential enters through Aij = |~ dr ■ A, 
which is a quenched random variable on each link in 
the lattice, with a uniform probability distribution in 
[0, 271"), corresponding to fixed random magnetic induc- 
tion B = V x A through the plaquettes of the lattice. 

The partition function is Z = J Q 27r d6j^ exp(— H/T), 

where T is the temperature. The simulation uses square 
lattices of size L x L with periodic boundary conditions 
(PBC) in both the x and y-directions. 

We also consider the GG model with different bound- 
ary conditions, where the system repeats periodically, 
but with fluctuating phase twist variables , A y added 
to the respective phase differences across the boundary 
of the system. This allows for phase twist fluctuations 
across the system, and enables calculation of the resis- 
tance of the model, as is discussed below. In this case 
the Hamiltonian is 

H = -J2 cos (^ - 0j - Aij - S Xi , L A x - S Vz xA y ) (2) 

(i,3) 

where A x , A y are integrated over (— oo, oo) in the parti- 
tion function. In a static simulation, the partition func- 
tion is invariant if the the phase twist at the boundary is 
spread out in the bulk as a constant extra phase difference 
A/L on each link. However, dynamically these both pre- 
scriptions are different. We define the phase twist fluctu- 
ations located at the boundary in order to maintain local 
dynamics. 

In the Monte Carlo simulations the trial moves are at- 
tempts to assign the phase, at a randomly chosen site, a 
random value from a uniform distribution in [0, 27r). The 
moves are accepted with probability 1/(1 + expAE/T), 
where AE is the energy difference for the trial move. One 
sweep is defined as on average one attempt to update 
each phase. In the simulations of static quantities we 
use a temperature exchange methodiS*^, that consider- 
ably enhances the convergence of the simulations at low 
temperature and big system sizes. We use ten sweeps 
of local MC moves, followed by one exchange sweep. 
In the exchange sweep the trial moves are attempts to 
exchange configurations between neighboring tempera- 
tures, where up to 39 closely spaced temperatures, with 
T n = T„_i(l — 1/L), are simulated in parallel, sharing the 
same disorder realization. The exchange moves are ac- 
cepted with probability exp[-(/3„ - (3 n+1 )(E n+1 - E n )], 
where (3 = 1/T. To reach equilibrium we discard up 
to 6 ■ 10 4 MC sweeps, followed by equally many sweeps 
where data is collected. The results are averaged over up 
to 10 4 realizations of the random vector potential A^. 
To avoid systematic errors in the calculation of squares 
of expectation values, two replicas of the system with 
the same disorder are used. We use various tests to ver- 
ify that sufficient equilibration is achieved, which will be 
presented in the results section below. In this paper we 
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denote thermal averages by (•••), and disorder averages 
by [.-.]. 

In the simulation of dynamic quantities the exchange 
MC method does not apply directly, since we must use 
a method that respects the local nature of the dynam- 
ics of the system. Instead we use the usual Monte Carlo 
scheme with local updates of the phase. However, also 
for the dynamics runs we take advantage of the exchange 
algorithm in two ways. First, we use equilibrium configu- 
rations from an exchange run as initial configurations for 
the dynamics runs. Second, during the dynamics simu- 
lation we also recalculate the static quantities, and com- 
pare to the same quantities from the exchange runs. This 
gives an additional test of convergence for the dynamics 
simulation. We define one time step to be one MC sweep 
through the system with local updates. In the simula- 
tion with fluctuating twists, one time step corresponds 
to one sweep with local phase updates, followed by one 
update of A x and of A y . The phase twist updates are at- 
tempts to add a uniformly distributed random number in 
(— 7r,7r). As initial configurations in the dynamics runs, 
we use configurations generated in an exchange run, fol- 
lowed by up to 4 • 10 5 sweeps of local MC moves where 
data is taken. Occasional runs with even more sweeps 
were also used to test convergence. Disorder averages are 
formed over up to 10 4 realizations of the disorder. We 
stress that there are many dynamic universality classes, 
and here we only consider that of local relaxation of the 
phase, which has proven useful in previous studies of the 
superconducting phase transitio n 13 ! 23 . 



III. CALCULATED QUANTITIES AND FINITE 
SIZE SCALING RELATIONS 

The basic quantity in our analysis of thermodynamic 
properties of the glass transition is the root mean square 
(RMS) currenlj 4 *^. The current in the x-direction is de- 
fined as the derivative of the free energy F with respect 
to a phase twist A x on the boundary in the x-direction, 
for A x = 0: 



dF 



dA x 



A.=0 
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(sin( 



'j+x 



4«+*)> ( 3 ) 



where j + x denotes the nearest neighbor (NN) site next 
to j in the ^-direction. Here the phase twist at the 
boundary has been absorbed into the bulk by adding a 
uniform phase twist A = A x /L on each link in the in- 
direction, which gives the 1/L factor in the last equation. 
The RMS current is defined as 



(4) 



where a, b denote two replicas of the system sharing iden- 
tical disorder. 

In order to use the scaling properties of MC data to 
discriminate between the possibilities that T c = and 
T c > 0, we must formulate suitable scaling assumptions 



for each cas o 3 i 13 . At the glass transition the correlation 
length is assumed to diverge as £ ~ \T — T c \~ v if T c > 0, 
and as £ ~ T~ v if T c = 0. In the case T c > 0, we assume 
that close to T c the RMS current obeys the finite size 
scaling relation^ 



I Ims (T,L)=I(L i /"(T-T c )) 



(5) 



where I is a universal scaling function. This scaling form 
implies that the RMS current is independent of system 
size L at T = T Cl which gives a convenient criterion for 
locating T c . In contrast, if T c — 0, the scaling relation 
becomes modified, since for T c — we have £ ~ T~ v = 
V , and hence ~ which must be included in the 

scaling relation. This gives (31 — d(3F/dA ~ £°. The 
scaling ansatz for the RMS current valid for the case 
T c = becomes 



I Tms (T,L) = L- 1 / v I(L 1 / v T) 



(6) 



Below we will compare both these scaling forms to our 
data, and find that the scaling form with T c > is not 
fulfilled by our MC data, while the scaling form with 
T c = works well. 

Next we discuss the quantities that we calculate in 
the dynamic simulation, and their respective scaling re- 
lations. The aim of our analysis is to obtain scaling rela- 
tions from a dynamic scaling ansatz that assumes T c = 0, 
and here we will only formulate the scaling relations for 
this case. The scaling relations will first be formulated 
in terms of the correlation time, r, whose scaling will in 
turn be discussed below. 

The resistance can be calculated by using fluctuating 
twist boundary conditions instead of PBC. It is obtained 
from the Kubo formula for the voltage- voltage correlation 
function. We define a time dependent resistance function 
as 



R(t,T,L) = ± £ At'[(V(t')V(0))} 



(7) 



where the voltage at MC timestep t is given by the 
Josephson relation V(t) = ^r^ 2 - (we henceforth assume 
units such that h/2e = 1). Here dA x /dt is the rate of 
change of the phase shift across the system in the x- 
direction, which gives the resistance in the ^-direction. 
The usual Kubo formula 2 ^ for the resistance is obtained 
from this function in the limit when the summation over 
MC timesteps t — > 00. In practice this limit means that 
the summation time has to be long enough that the re- 
sistance is independent of t. However, since this time is 
going to be very long in the simulation, it is quite useful 
to consider R as a function of t and study the scaling 
properties as function of t. Since the voltage scales as 
V ~ 1/r, we make the finite size scaling ansatz 



R(t, T, L) 



±.R(L^T,t/r) 



(8) 



where R is a scaling function. 
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Another useful quantity is the dynamic correlation 
function for the real part of the superconducting order 
parameter, 9' — 5Z,-cos0j, which is gauge dependent 
and hence not directly observable in a superconductor, 
but gives a convenient measure of the dynamic correla- 
tions of the local order in the system. We define the time 
dependent summed autocorrelation function 



G(t,T,L)= ]T At' 



(9) 



For t — > oo, G measures the correlation time r of the 
fluctuations in and we have 



G{t, T, L) = tG{L x I v T, t/r) 



(10) 



Note that all the (unknown) dynamic scaling functions 
contain two arguments, which complicates their analysis. 
We will make use of two possible ways to obtain a func- 
tion of one variable only, and this function can then be 
conveniently analyzed. One way is to obtain data for 
t t, and the second way is to fix the first argument to 
a constant value by using different temperatures for each 
system size so that L X I V T = const 1 ^, which then permits 
scaling analysis of the entire t/r dependence. 

In order to quantitatively fit the dynamic scaling for- 
mulas to MC data, we need an ansatz for the form of 
t. This is related to the divergence of the free energy 
barriers at the glass transition. The thermodynamic free 
energy barrier exponent 8 is defined from the interfacial 
free energy RMS difference, Ul, for changing the bound- 
ary condition in one direction from PBC to anti-PBC. 
For a system of size L we have Ul = KL , where K is a 
constant. Now, a change in free energy from a change in 
the boundary conditions must scale in the same way as 
the current in Eqs. © and (JBJ, which gives 8 = —1/v. 
Thus 8 can be obtained from analysis of / at finite tem- 
peratures. 

A dynamic barrier exponent ip can be defined by Vl = 
KL^, which is the free energy cost to move a vortex a 
distance L using only local moves. In general Vl and 
Ul do not need to coincide, and the dynamic barrier for 
vortex motion across the system can be very big. We 
assume that the correlation time scales in an activated 
form given by 



t ~ exp(VyT) ~ exptC/T 1 ^) 



(11) 



where C is a constant. If ip = this reduces to a simple 
Arrhenius form, t ~ expC/T. We also consider a few 
other possibilities. The usual finite temperature form of 
the critical slowing down is r ~ L z . Another possible 
form is a logarithmic barrier, r ~ exp K (In L)/T. 



IV. RESULTS FOR STATIC QUANTITIES 

We first study the approach to equilibration of the ex- 
change MC simulation. This is necessary in order to 




FIG. 1: Test of equilibration of MC data for the RMS current 
at T = 0.1. Here t sweeps are taken without any measure- 
ments, followed by another t sweeps where data is collected. 
The saturation time in the plot gives an estimate of the nec- 
essary equilibration times for the different system sizes. 



avoid systematic errors due to insufficient equilibration. 
We consider the RMS current, I lms , defined in Eq. 
Figure^shows some of our data for the RMS current as a 
function of simulation warmup time. The data points in 
the figure are calculated by first discarding t MC sweeps, 
and then forming averages over equally many sweeps. In 
Fig.^we see that the equilibration time, where the curves 
become independent of t, becomes very long already for 
moderate system sizes L, despite the considerable con- 
vergence speed-up due to the exchange method. 

Next we consider finite size scaling properties of the 
RMS current. Figure H shows MC data for the RMS cur- 
rent vs. temperature. According to Eq. data curves 
for J rms should become system size independent and thus 
intersect at T = T c , if T c > 0. The absence of any such 
intersection in the figure indicates that T c is much smaller 
than the temperatures where we have data, i.e., T c < 0.1. 
Henceforth we will therefore assume that T c = 0. 

To scale MC data for J rms assuming that T c = 0, Eq. 
© shows that a plot of L 1 /"/,™ vs. L X I V T for differ- 
ent T, L should collapse onto a system size independent 
curve, for the correct value of \jv. Scaling plots are 
shown in Fig. |3| for the choices (a) \jv — 0.39 and (b) 
1/u = 0.36. In (a) we see that 1/v = 0.39 gives the best 
overall data collapse, but the fit is somewhat poor for the 
smallest values of L 1 / i/ T, which for T c = is expected 
to be the most significant part of the data collapse. We 
therefore drop data points for large L 1 / !y T and only fit 
the value of v to the data points with smallest L X I V T . 
Fitting the data points at the two lowest temperatures 
for sizes L = 10, 12, 16 to a straight line, and minimizing 
the RMS fit error, gives v = 0.36 ± 0.03. The error bar 
is obtained by the bootstrap method^. The resulting 
data collapse is shown in Fig. 0(b), and we see that, for 
small L 1 /' y T, a better fit is obtained than in (a). The es- 
timate v = 0.36 ± 0.03 is insensitive within the error bar 
to varying the precise details of the fit, such as including 
the three lowest temperatures for each L. 
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FIG. 2: MC data for the RMS current vs. temperature for 
different system sizes L. Finite size scaling arguments valid 
for T c > (see text) show that data curves for different system 
sizes become independent of L at T — T c , and hence intersect. 
The absence of any intersection in the plot implies that T c <C 
0.1. The data points connected by lines are produced by the 
temperature exchange method. The circles mark data points 
obtained without the exchange method, to allow calculation of 
dynamic quantities, and are in good agreement with the data 
points obtained from the exchange method, indicating that 
ensemble averages and dynamic averages tend to coincide. 
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FIG. 4: Examples from our simulation of the time evolution 
of MC data for the phase twist variable for a few different 
disorder realizations. The characteristic time scale for big 
fluctuations in the twist variable can be > 10 sweeps. We 
identify such fluctuations as phase slips corresponding to vor- 
tex motion across the sample. 



V. RESULTS FOR DYNAMIC QUANTITIES 
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FIG. 3: Finite size scaling plot of the MC data in Fig. [5] 
using Eq. © with T c = 0. In (a) l/u = 0.39, which gives a 
data collapse except for the lowest temperatures and smallest 
system sizes. In (b) l/u — 0.36, which gives a better data 
collapse at the smallest temperatures. 



The dynamics simulations contain only local MC 
moves, and the averages computed in the dynamics runs 
correspond to time integration of the observables. The 
statics simulation contains both local MC moves and non- 
local temperature exchange moves that take big steps in 
configuration space, which corresponds to ensemble av- 
eraging. If ergodicity is broken these different averages 
may disagreed. To investigate this issue we plot the RMS 
current calculated in a dynamics run in Fig. [21 together 
with the corresponding result from the static simulation. 
The figure shows that the results from static and dynamic 
runs coincide within error bars. 

We will make use of the following input from the statics 
simulations. The static scaling results gives the value of 
the correlation length exponent, which is needed also in 
the dynamic scaling analysis. We will assume the value 
l/u = 0.39, but make occasional comparison with the 
corresponding results for l/u = 0.36. Furthermore, the 
static scaling results also give an indication of the range 
of system sizes and temperatures where we can expect 
dynamic scaling of the data to hold. From Fig. [3] we esti- 
mate that dynamic scaling can be expected when L 1 / !y T 
is smaller than about 0.8. 

We expect that the longest emergent time scale comes 
from phase slip events corresponding to vortex motion 
across the sample. An illustration of these long corre- 
lation times is seen in Fig. 0] which shows phase slip 
events in a dynamic simulation where fluctuations in the 
boundary conditions are included, for a few different re- 
alizations of the quenched disorder. Phase slip events are 
often as far as about 10 5 sweeps apart. 

Figure shows the time dependent resistance func- 
tion given by Eq. Q for L 39 T = 0.8. The data curves 
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FIG. 5: MC data for the resistance function R(t,T,L) vs. 
time t for different system sizes for A = L 1/0 ' 39 T = 0.8. The 
resistance curves have different characteristic regimes, one for 
t > t where R becomes independent of t, and one power-law 
regime for shorter times. 
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FIG. 6: MC data for RTt vs. t/r for different temperatures 
and system sizes with A = L ' 39 T = 0.4 and A = 0.8. For 
each value of T, L the correlation time r has been tuned until 
an optimal data collapse is obtained in the large t limit. 



show unexpectedly complicated behavior, with a power- 
law regime at short times, and a time-independent limit 
at long times, with various crossover times in between. 
We will now analyze the scaling behavior of these differ- 
ent time scales. 

We first consider the long time limit, t 3> r, of the data 
curves in Fig.0 which gives the resistance R. According 
to Eq. © , the correct value of r for each system size and 
temperature, such that L X I V T = A is a constant, should 
make data curves for RTt plotted vs. t/r collapse onto a 
common curve for f>r. Thus we adjust r for each data 
curve to the value that gives the best collapse. The result 
is shown in Fig. HJ1 We obtain a data collapse to a com- 
mon curve at large t, for each set of data curves belonging 
to the same value for A — L X I V T . At shorter times the 
collapse breaks down, indicating a different short time 
behavior. By this procedure we determine an estimate 
of the correlation time r(T, L), up to an unknown mul- 
tiplicative constant that depends on A — L X I V T . The 
scaling result shown in this figure demonstrates that, up 
to the level of statistical accuracy of the simulation data, 
the resistance obeys a dynamic scaling relation which as- 
sumes T c = 0. This result is independent of any assump- 
tion about the functional form of t(X, T). 

To proceed further, we need to fit the values for r 
and we then need an assumption of a functional form 
for r. All the functional forms r ~ exp K (In L)/T, 
t ~ cxpKL^/T, and r ~ L z approximately fit the 
data, and our present data is therefore not sufficient 
to clearly distinguish between these functional forms. 
Larger system sizes and lower temperatures are neces- 
sary in order to discriminate between them. The form 
r ~ exp_ftT(lnL)/T gives a good fit to the data, with 
K « 0.56 independent of T and L, shown in Fig. [7| 
However, this should only be considered as one possible 
parameterization of the correlation time. 

Next we consider the scaling behavior of R(t, T, L) for 



■ 


A = 0.4 






1 






▲ 


A = 0.8 
























































..it'' 


































1 







5 10 

T 1 In L 



FIG. 7: MC data for the correlation time t, determined in 
Fig.|H] vs. (In Li) /T for different temperatures and system sizes 
with A = L°' 39 T = 0.4 and A = 0.8. The dashed lines corre- 
spond to t ~ exp K (In L)/T, with K = 0.56 independent of 
T,L. 



intermediate times. Surprisingly the scaling found for 
long times in Fig. is not obeyed for intermediate times 
where the data collapse breaks up. We therefore attempt 
another scaling relation to apply in this regime. Accord- 
ing to Eq. |JSJ|, we again set L X I V T = A to a constant 
value, and plot RTt vs. t/r. We can then obtain an ap- 
proximate data collapse in the intermediate power-law 
regime in the figure for r ~ L z with z — 2.7 ± 0.2, 
where the error bar is estimated from the interval out- 
side of which a significantly worse fit is obtained. Wc 
can also fit the same data with other functional forms 
for r, but we get the best fit with r ~ L z . These re- 
sults show that the resistance has two separate dynamic 
scaling regimes, characterized by two separate diverging 
correlation times. It is natural to associate the longest of 
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FIG. 8: Finite size scaling collapse of MC data for the time 
dependent resistance function R(t,T,L) with system sizes L 
and temperatures T related by A = L 1/u T = 0.8. The data 
collapse is obtained for r ~ L z with z = 2.75. Data for other 
values of A scale with nearly the same value of z. 
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FIG. 9: Finite size scaling data collapse of MC data for the 
dynamic correlation function G(t,T,L) defined in Eq. 0. 



these correlation times with the phase slip processes that 
make the resistance finite at any finite temperature, since 
the correlation time roughly matches the time between 
phase slips in Fig. 0] We will now study these different 
dynamic scaling regimes in more detail. 

To reach further insight into the scaling properties at 
intermediate times we consider the order parameter auto- 
correlation function G(t, T, L), defined in Eq. @. Figure 
shows MC data for G(t, T, L) scaled according to Eq. 
JHJ with t ~ L z . The data in the figure has L X > V T = 0.8, 
and the data collapse is obtained for z — 2.74. The data 
in the figure is for the case of fluctuating phase twists 
at the boundaries. We also calculated the same quan- 
tity with periodic boundary conditions, and obtained an 
equally good fit with the same exponent. Approximately 
the same exponent also gives a collapse for other values 
of A, but a slight drift in z is noticeable, just as the data 
collapse in Fig. [31 for the RMS current is somewhat de- 
pendent on A. Including this uncertainty we estimate 
z = 2.7 ± 0.2. We conclude that the dynamics measured 
by G is independent of the phase twist fluctuations. The 
dynamic scaling relation r ~ L 2 ' 7 , instead comes from 
local dynamic fluctuations other than the global phase 
slips. This means that we can interpret the intermediate 
resistance scaling regime in Fig. as coming from local 
fluctuations rather than global phase slips. 



VI. DISCUSSION AND CONCLUSIONS 

This paper considers finite size scaling properties of 
Monte Carlo data for the glass transition in the two- 
dimensional gauge glass model. Equilibration of our 
simulations at low temperatures is achieved by a tem- 
perature exchange simulation method. Up to the level 
of statistical uncertainty, system sizes and temperatures 
reached in the simulations, our finite size scaling results 



for both static and dynamic properties show scaling prop- 
erties that are clearly consistent with a glass transition 
occurring only at zero temperature, and thus indicates 
the absence of a superconducting glass state at any finite 
temperature in 2D. 

Our results for static scaling properties reproduces, 
and adds some details to, some of the results of Refs. 
EillOl in addition to providing a useful starting point for 
our dynamics simulations. From the scaling of the RMS 
current we find the correlation length exponent given by 
l/V » 0.39, obtained by optimizing the quality of a fi- 
nite size scaling data collapse over a finite temperature 
interval close to T = 0. This value agrees with that of 
Refs. HlH If we instead optimize the fit at the lowest 
temperatures where we have data, where scaling is also 
expected to work best, we obtain 1/v = 0.36 ± 0.03. It 
is interesting to note that this is in very good agreement 
with the result for the ground state domain wall exponent 
9 = — 1/v in Ref. |8j. To obtain an even more accurate 
estimate of data for lower temperatures and bigger 
system sizes is necessary. To converge such data would 
require a considerable increase in simulation efforts. 

For the dynamic properties, we consider the linear re- 
sistance computed from the diffusion of the phase slip 
across the boundary of the system. Since the correlation 
times grow very long at low temperature, it is not pos- 
sible to converge the simulation at as low temperature 
as for the static data where the temperature exchange 
algorithm helps considerably, which makes it more dif- 
ficult to study dynamic scaling than static scaling. A 
further complication is that our data for the linear resis- 
tance shows more than one characteristic time scale. We 
identify the longest correlation time as the characteristic 
time scale for phase slip fluctuations across the sample. 
Unfortunately it is not possible to determine the corre- 
lation time uniquely from the data, since the interval of 
sizes and temperatures where we have data is too lim- 
ited. All the forms r ~ expKL^/T, r ~ exp K In L/T, 
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and t ~ L z can be made to approximately fit the data. 
To discriminate between the possible forms, data at lower 
T and larger L is needed, which is beyond the present cal- 
culation. However, a logarithmically growing phase slip 
barrier, with r ~ exp K In L/T, gives a good fit to the 
data with K ks 0.56. A qualitatively similar result was 
obtained in a zero temperature calculation of the phase 
slip barrier energy in Ref. Qg, Since R ~ exp —KL^/T, 
of which the logarithmic barrier is a special case with 
ip = 0, a superconducting glass with zero resistance is 
obtained only for T = 0. In other words the dynamic 
glass transition occurs at zero temperature. 

Surprisingly, data for the time dependent resistance 
function R(t, T, L) does not scale with the longest cor- 
relation time for intermediate times. Instead, in this 
regime scaling is obtained by assuming a new correla- 
tion time of the form r ~ L z with z = 2.7 ± 0.2. Hence 
the scaling behavior of our data implies two different di- 
verging correlation times with quite different properties 
at the zero temperature glass transition. A further in- 
dication of the smaller correlation time is given by the 
scaling of the order parameter autocorrelation function, 
which scales over the entire time interval where we have 
data, including the longest times, quite in contrast to 
the resistance. The interpretation of this shorter correla- 
tion time is that it gives the characteristic time scale for 
fluctuations other than the global phase slip fluctuations 
that lead to a finite resistance. Similar results showing 
different dynamics for local fluctuations and phase twist 
fluctuations have been obtained previously at the finite 
temperature transition in XY models without disorder—. 
We also note that the scaling relation r ~ L z can be 
written r ~ exp z In L ~ expzL ~ exp C/T 1+ ^ u , which 



means that tp = 6 = —1/v. Hence the relation r <~ L z 
suggests that the static and dynamic free energy barrier 
exponents are the same for the local fluctuations of the 
order parameter. 

The two different diverging correlation times sug- 
gested here may explain why signatures of a glass phase 
with zero resistance at finite temperatures can easily 
be obtained in simulations. For the shorter correlation 
time we have i/> = 6 = which gives R(T) ~ 

exp— C/T + v v = const, for the resistance at small tem- 
peratures, i.e. the resistance stays at the T — value 
R = for small T. This shows that in a regime of inter- 
mediate time scales, zero resistance can be obtained at 
finite temperature. However, at much longer time scales 
that may be hard to reach in a simulation, phase slip fluc- 
tuation set in, and qualitatively different behavior follows 
with a superconducting glass transition only at zero tem- 
perature. 

In summary, from finite size scaling of a time depen- 
dent resistance function and an order parameter autocor- 
relation function, we constructed two different diverging 
correlation times from the T c = scaling behavior of 
Monte Carlo data. For further work, it would be of in- 
terest to further clarify the precise details of dynamic 
scaling at the zero temperature glass transition, and to 
include quantum fluctuations in the model. It would fur- 
thermore be interesting to search for the proposed novel 
dynamic scaling behavior in experiments. 
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